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Abstract 

We calculate the ratios i?^^^ = T{P eUe['y])/T{P fii^M {P = vr,/C) in 

Chiral Perturbation Theory to order e^p^. We complement the one- and two-loop 
effective theory results with a matching calculation of the local counterterm, per- 
formed within the large- A'^c expansion. We find = (1.2352 ± 0.0001) x 10"^ and 

R^^^ = (2.477 lb 0.001) X 10~^, with uncertainty induced by the matching procedure 
and chiral power counting. Given the sensitivity of upcoming new measurements, 
our results provide a clean baseline to detect or constrain effects from weak-scale 
new physics in these rare decays. As a by-product, we also update the theoretical 
analysis of the individual tt{K) — > ii>£ modes. 
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1 Introduction 



The ratio R\il = r(P ez/e[7])/r(P — > fJ^i^fj.i'j]) {P = tt,K) of leptonic decay rates 
of light pseudoscalar mesons is hehcity-suppressed in the Standard Model (SM), due to 
the V — A structure of charged current couplings. It is therefore a sensitive probe of 
all SM extensions that induce pseudoscalar currents and non-universal corrections to the 
lepton couplings p. Recently, attention to these process has been payed in the context 
of the Minimal Supersymmetric Standard Model, with [2] and without [3] lepton-flavor- 
violating effects. In general, effects from weak-scale new physics are expected in the range 
(ARe/fj) /Re/fj. ~ 10"^ — 10~^ and there is a reahstic chance to detect or constrain them be- 
cause of the following circumstances, (i) First, ongoing experimental searches plan to reach 
a fractional uncertainty of (Ai?^/J^)/i?^/^ ^ 5 x 10"^ [4] and {Ar[^^)/r[^^ ^ 3 x lO'^ [5], 
which represent respectively a factor of 5 and 10 improvement over current errors [6]. 
(ii) At the same time, the SM theoretical uncertainty can be pushed below this level, 
since to a first approximation the strong interaction dynamics cancels out in the ratio 
Re/^ and hadronic structure dependence appears only through electroweak corrections. 

Indeed, the most recent theoretical predictions read -R^^^ = (1.2352 ± 0.0005) x 10~^ [7], 

^ (1.2354 ± 0.0002) x 10"^ [H], and R^J^ = (2.472 ± 0.001) x 10^^ [8]. In Ref. [7j a 
general parameterization of the hadronic effects is given, with an estimate of the leading 
model-independent contributions based on current algebra |9]. The dominant hadronic 
uncertainty is roughly estimated via dimensional analysis. In Ref. [8], on the other hand, 
the hadronic component is calculated by modeling the low- and intermediate-momentum 
region of the loops involving virtual photons. 

The primary goal of this investigation is to improve the current status of the hadronic 
structure dependent effects. To this end, we have analyzed Re/fj, within Chiral Perturbation 
Theory (ChPT) [10], the low-energy effective field theory (EFT) of QCD. The key feature 
of this framework is that it provides a controlled expansion of the amplitudes in terms 
of the masses of pseudoscalar mesons and charged leptons {p ~ m^ ^^^^/A^, with A^ ~ 
AttFj^ ~ 1.2 GeV), and the electromagnetic coupling (e). Electromagnetic corrections to 
(semi)-leptonic decays of K and vr have been worked out to O(e^p^) pTl[T2] . but had never 
been pushed to O(e^p^), as required for Re/fi in order to match the experimental accuracy. 
In this work we report full details of our analysis of Re/fj, to O(e^p^), while a summary 
of the results is presented elsewhere [13]. To the order we work in ChPT, R^/^ features 
both model independent double chiral logarithms (previously neglected) and an a priori 
unknown low-energy coupling (LEG). By including the finite loop effects and estimating 
the LEG via a matching calculation in large- A^^c QGD, we thus provide the first complete 
result of Re/fj. to O(e^p^) in the EFT power counting. Most importantly, the matching 
calculation allows us to further reduce the theoretical uncertainty and put it on more solid 
ground. 

Our presentation is organized as follows. In Section [2] we introduce the basic definitions 
and outline the strategy to calculate Re/fj. to O(e^p^). In Section [3] we shortly review 
the basic GhPT formalism and the needed effective lagrangians. The loop calculation is 
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described in Section |4] and in Appendix [XJ and the results are reported in Section [51 We 
then report the matching calculation of the effective coupling in Section [6l with technical 
details in Appendix [Bl We present the contribution from real photon emission in Section [3, 
while in Section [8] we give our final analytical and numerical results for -Rg/^^'' and discuss 
them. Section[9]is devoted to updating the theoretical expression for the individual 7r{K) — >• 
iui rates. Finally, Section [10] contains our concluding remarks. Since we are reporting here 
the first ChPT calculation to order e^p^, we give several details and intermediate steps of 
our analysis, both throughout the text and in the Appendixes. 



2 R^Y^^ in ChPT: overview 

To avoid excessive notational clutter, throughout this paper we illustrate the main argu- 
ments in the case of vr — s> £z/ decays and subsequently report any significant changes that 
occur for K decays. We consider the ratio 

dW ^ r(7r+^e+;/e(7)) 

to order e^p^ in Chiral Perturbation Theory (ChPT). Within ChPT the invariant ampli- 
tudes can be expanded in powers of the external masses and momenta (of both pseu- 
doscalar mesons and leptons) and powers of the electromagnetic coupling. To leading order 
in the chiral expansion one finds 

Tf = -i2GfV:^F m, ul{Pu) v{pi) . (2) 

F can be identified to lowest order with (and Fk^F^)- Setting e = 0, to a given 
order (p^*^) in the purely "strong" chiral expansion, the amplitude reads as above, with 
the replacement F — ^ Fn'^^\ Fi^"^ being the pion decay constant to order p^". When 
considering the ratio of electron-to-muon decay rates the pion decay constant drops and 
one obtains the well known expression: 

Non-trivial corrections to Eq. [3] arise only when e 7^ 0, i.e. to order e^p^" in ChPT. 

Lorentz invariance implies that higher order contributions are proportional to the lowest 
order amplitude, and this allows one to write to O(e^p^) 

r(7r ^ tuYl]) = r(°)(7r ^ iu) X [1 + 2 Re [rf^" + rf + 5f^" + 5f^"] , (4) 



^ Intermediate steps in our analysis depend on the definition of the invariant amphtude Ti {£ = fi,e), 
for which we use out{i~^{pi)i^e{Piy)\'^^{p))in = (27r)'*(5('*' {p - Pe - Pv) iTt. 
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where 



and 



r(o)(7r^£i/) 



Air 



■ rrtTT rrii, 



m 



mi 



2\ 2 



TP 

2p2n _ r(7r — » fz/7)|e2p2r, 

0* — 



r(0)(7r ^ £z/) 



(5) 

(6) 
(7) 



are respectively the corrections induced by virtual and real photon effects, whose sum is 
free of infrared divergences. Taking the ratio of electron and muon decay rates one obtains: 



ejii 



R 



e/fi 



^2p4 -r 



2Re(r:V"-<^") + (5:V"_5e 



(8) 



(9) 



The main feature emerging from Eq. M is that only those diagrams that depend in a 
non-trivial way on the lepton mass contribute to -Re/^t- The diagrams leading to me- 

g2 2n 

independent will drop when taking the difference of electron and muon amplitudes. 

This observation greatly reduces the number of diagrams to be calculated in the effective 
theory. All the considerations presented in this section trivially extend to the case of 
leptonic decays of charged kaons {K — > iu). 



3 Electromagnetic corrections to (semi)-leptonic pro- 
cesses at low energy 

The appropriate theoretical framework for the analysis of electromagnetic effects in semilep- 
tonic kaon decays is a low-energy effective quantum field theory where the asymptotic states 
consist of the pseudoscalar octet, the photon and the light leptons [TT]. The corresponding 
lowest-order effective Lagrangian is given by 

+ ^[l(z^ + e4-m,)£ + I^z^z/,L]. (10) 

e 

F denotes the pion decay constant in the chiral limit and in the absence of electroweak 
interactions. The low energy constant Z ~ 0.8 can be determined by mass splitting of 
charged and neutral pions. The symbol ( ) denotes the trace in three-dimensional flavour 
space, and 

= i[u\d^ - ir^)u - u{d^ - il^)v)] , (11) 
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with the Goldstone modes collected in the field u: 



u 



exp 



2$ 



V2 



TT 

K- 



1 



TT 



"76^8 



'12) 



The photon field and the leptons i^^ = e, /i) are contained in ( fTTl) by adding 
appropriate terms to the usual external vector and axial-vector sources v 



I, 



wf > 



(13) 



The 3x3 matrices Qi\i Qj^ are spurion fields. At the end, one identifies Qi\ with the 
quark charge matrix 

" 2/3 " 
g°'^= -1/3 , (14) 
-1/3 



whereas the weak spurion is taken at 

QI = -2V2Gf 



Vud 






(15) 



where Gp is the Fermi coupling constant and Vud, Kts are Cabibbo-Kobayashi-Maskawa 
matrix elements. For the construction of the effective Lagrangian it is also convenient to 
define 

Q^--" := nQ^'"«^ QT 

Explicit chiral symmetry breaking is included in x+ = 
to the quark mass matrix: 



tn"^„ (16) 
ux^u where x is proportional 



X 



2Bn 










rrid 
m. 



(17) 



and the factor Bq is related to the quark condensate in the chiral limit by (0|gg|0) = —F^Bq. 
The local action at next-to-leading order involves the sum of three terms, Cp* + -^^2^2 + 

'^e*2pV The first one, Cp4 includes the well-known Gasser-Leutwyler Lagrangian [13] in the 
presence of the generalized external sources introduced in f[T5]) . as well as a term from the 
Wess-Zumino-Witten functional that incorporates the effect of chiral anomalies [TS]. Here 
we quote only the operators relevant to our analysis: 



£p4 



487r^ 
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(a) (b) 
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Figure 1: Diagrams contributing to Re/^ to order e^p^. Dashed lines indicate pseudoscalar 
mesons, solid lines leptons, and wavy lines photons. 
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The second term, C^Jipi, encodes the interaction of ultraviolet (UV) virtual photons with 
hadronic degrees of freedom [16], [ITJ [H]. It contributes to the individual P ^ eu and 

2 2 

P — > ^p, but leads to an m^-independent r| ^ so that it has no effect on Re/^- The 
same argument applies to £'T*2 , which involves leptonic bilinears. Similarly, when inserted 
in one-loop purely mesonic graphs, these effective operators contribute to P ^ eu and 
P ^ fiu to order e^p^, but their contribution cancels in Re/fi- Therefore, there is no need 
to report the full expression of these effective lagrangians here. 

Finally, we shall see that a counterterm of O(e^p^) is needed in order to make Re/fj. finite 
to O(e^p^). While we have not constructed the most general >C|f2p4, on the basis of power 
counting we can conclude that the same combination of operators (and LECs) contributes 
to both R[^^ and -Rg/^- This fact is also explicitly borne out in the matching calculation 
that we perform in Section [6l 



4 Virtual-photon corrections: analysis 

We work in Feynman gauge, use dimensional regularization to deal with ultraviolet (UV) 
divergences and an infinitesimal photon mass to deal with infrared (IR) divergences. We 
report the diagrams contributing to Re/^ to 0{e^p^) and 0{e^p^) in Fig. [T]and Figs. [2J131 
respectively. At the order we work, we need the charged lepton and pseudoscalar meson 

2 

wavefunction renormalizations to one-loop accuracy. We denote them hy = I + 
(charged lepton) and = 1 + + (pseudoscalar meson). 
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Figure 2: IPI diagrams contributing to R^/^ to order e^p^. Shaded squares indicate 
vertices from the O(p^) effective lagrangian. 

To order e^p^ one has to consider only two IPI diagrams and the effect of charged lepton 

2 2 

wave-function renormahzation (see Fig. [1]). The resulting amplitude T'^ ^ [11] coincides 
with the point-like approximation of Ref. [19]. Since this is well known, we do not dwell 
further on it, but we will simply report the result in the next section. The situation is 
more interesting to next-to-leading order. 

4.1 Organizing the 0(e^j9^) diagrams 

To O(e^p^) one has to consider (i) two- loop graphs with vertices from the lowest order 
effective lagrangian and (ii) one-loop graphs with one insertion from the NLO lagrangian 
£p4 (we denote the latter vertices with shaded squares); (iii) a tree level diagram with 
insertion of a local operator of O(e^p^). In Fig.[2]we report all relevant IPI topologies: each 
O(p^) vertex receives contributions from several 0{p^) operators and all allowed mesons 
run in the internal loops. External leg corrections are depicted in Fig. [31 

The self-energy insertion on the internal mesonic leg (class {IV) in Fig. [2l) is handled by 
observing that to O(p^) the self-energy reads S(p^) = A + Bp^ (with A and B momentum- 
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Figure 3: External leg corrections to Re/^ to order e^p^. 



independent) and therefore 




d 



(20) 



2 ' 

'0 



where Z represents the on-shell wave-function renormalization, mo is the 0{p'^) mass and 
m is the physical 0{p^) mass. With this result at hand, by re-grouping the diagrams of 
class {IV) and external leg corrections with those of classes (/), (//), and {III), it is 
straightforward to show that the inclusion of virtual corrections to 0(e^j9^) amounts to: 

• using the physical O(p^) meson mass in the amplitude of O(e^p^); 

• calculating a set of "effective" one-loop diagrams with vertices given by appropri- 
ate off-shell form factors evaluated to 0{p^) in d-dimensions. These effective one- 
loop diagrams are shown in Fig. HI The shaded circles denote respectively: the 
d-dimensional 0{p^) iriu vertex (Fig. Sl^a) and (d), with off-shell pion and charged 
lepton in Fig. |4](a)); the d-dimensional O(p^) vr7r7 vertex with the photon and one 
pion off-shell (Fig. ID^b)); the d-dimensional O(p^) 7rii"~f vertex with the photon and 
charged lepton off-shell (Fig. 111(c)). 

Within this approach one starts the calculation of genuine two-loop diag rams at a stage 
where the one-loop sub-divergences (generating non-local singularities) have already been 
subtracted. As we shall see, another advantage is that the non-local 0{p^) vertices admit 
a simple dispersive parameterization that greatly simplifies the calculation. 

As seen from Fig. |H the virtual photon contributions can be divided into IPI and 
external leg corrections. For the external leg corrections we find: 




non-lP/ 



1 




(21) 
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where F^^'^ / F has to be evaluated in d-dimensions. The IPI contribution can be written 
as the sum of the mass-renormahzation in Tf ^ and a convolution: 



2Gj.V:,e'F 



ipi 



(22) 



where 



T, 



V-A 



V2F 



dx e^^^+^^J^ (0|T(J^^'^(x) (K - A,){y)\7r+{p)) 



(23) 



with V^{Afj_) = u'-ff^{'-f^)d and W = p — q. Lorentz invariance and Ward identities imply 
that T^i^~^ in turn can be decomposed as follows (see also [20]): 



(T^-^)'^^(p,g) = ^V,e'^•'-^q^pp + 



{2p — qY{p — qY 



I-IU 



F^' 



(4) 



2p ■ q — q"^ 

- A^{q- pg^'^ - p^q^) - {A^ - A^) {q'g^'' - q^q'') 
{2p- qY{p- qf qf'ip-qy 



2p ■ q — q^ 



{Fy'^iq') - 1) 



- As[q-p {q^p" - q^q") + (p^g" - p^'p")] 



(24) 



The form factors Vi, Ai depend in general on both q^ and W"^ = {p — and have to be 
evaluated to O(p^) in ChPT in d-dimensions El. The same applies to the pion form factor 
Fy^{q^) and decay constant F^. The convolution integral generates a term proportional to 

2 2 

^ I IP/ as well as terms induced by V^i, 2, and Fy"^ — 1. With obvious notation we can 
write 



Tie 



IPI 



+ 



/ 2 I 9 1 

IP/ V "Ip" "IpV « 



•9 

J- n 



(25) 



Combining Eqs. [25] and [21] we then obtain: 



(4) 



1 1 Tf^" 



m. 



TT U4 



m 



d 



dmi 



■T, 



e-p- 



(26) 



The effect of the the last two terms in Eq. [26] is taken into account by simply using the 

2 2 

physical pion mass and decay constant to 0{p^) in ^ . The remaining terms provide a 
genuine shift to the invariant amplitude. In order to calculate such a shift, we need to: 



^In this work we use the convention £0123 = +1 for the Levi Civita symbol. 
■^To 0{p*) the form factor A3 vanishes. 
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(b) 
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Figure 4: Effective one-loop diagrams contributing to Re/^ to order e^p^. The shaded 
circles represent the 0{p^) contribution to d- dimensional off-shell effective vertices. 

(i) Work out the relevant form factors Vi, Ai^2, Fy^ to O(p^) (one-loop) in d-dimensions. 

(ii) Insert them in the convolution representation of Eq. [22] and calculate the resulting 
integrals. 

In the following subsections we report the results of these steps. 
4.2 Form factors in d-dimensions 

We work in d dimension with d = 4 + 2w [2T1I22]. The relevant form factors to O(p^) read: 

Nc 



Vi = 
A, = 

A2 = 

F^^{q^) = l + 2H^^{q^) + HKK{q') 
The loop function Haa{q^) [E] reads 



247r2F2 

4 (Lg + Lio) 

2 



with 



A{m^) 

r^iq^] 
r^iq^) 



A{ml) d-2 2 
ml 8{d-l) ^ 3 ^ 



d'^k 



(f-^rnl - 



A{d-V 



r'^iq^) 



m 



2+2w 



{27rY ^^-711? + ie {Air] 



2+w 



rr-1 - 



—I 



d'^k 



{27tY -ml + ie] [{k - qY -ml + it] 



(4vr; 



2+w 



I dx [m^ — q^ x(l — x 
Jo 



(27) 

(28) 

(29) 
(30) 

(31) 

(32) 
(33) 

(34) 
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The function J°'°'{q^) admits a dispersive representation [211 122] in d- dimensions, which 
proves very useful in the evaluation of genuine two-loop contributions: 



[da] 



m 



2w 



[da] 



Uml a — q 

da r (I) 



(47r)2+"' r (I + w) \^ml 



a 



Ami 



2\ 1/2 



a 



(35) 
(36) 



4.3 "Effective" one- loop diagrams 



The IPI contributions Ty^ , T^^ , , Tp^ can be written as the convolution of a known kernel 
with the d-dimensional form factors Vi , 2, Fy'^ . It is simple to check that Ty^ , T^^ , , Tp^ 
are IR finite, so we set = 0. Upon inserting the O(p^) ChPT form factors into Eqs. [2l] 
and [22l we obtain: 



A(mj 



2{d-l)mj 
Aimj) 



4{d-l)mj 



[(4 + d)m', -{d- 2)ml] 
[(3ci2 -6d + 4:)m] - {d - 2fml] 



(37) 
(38) 



-2T; eH d (2a,, + Oi^i^) A(m|) 



+ 
+ 



26. 



8m: 



2/, 



(^)7r7r 



2 t-Wtttt 



2 Ai)KK 



d 



2 (/r - ID + (4^^ - ^^^) 



T7T7T 
-'3 



4m^ ( h 



KK 



(39) 



2TfeM(2a„ + a,,)!^i'^('"-'-'"'^('"'' 



ml — mj 



+ 



?72_ TTln t \ / 



m^ 6. 



m^ — m^ z 



TTTT / o j-{TT)mT . j(it)KK 



+ {ml + miy— (2Tr + ^ - SmlT^ - 4m^T/^^) 



(40) 



11 



In the above expressions we have used the definitions: 



1 [A{ml) d-2 , 2^ \ 
= F-- \ ml 8(^-1) + 3^V ^ ^ 

1 M(m2^) rf-2 2 \ 

*" = 4(d-l)P 

which come from the decomposition i/mm (9^) = Omm?^+^mm(Q'^— 4m,^) J™™(g^). Moreover, 
the building-block two loop integrals are defined as follows: 

j{£)aa 

-'1 = 



j{7T)aa 
-'1 



j[l)aa 
-'2 



^2 



{'K)aa 



jaa 
^3 



jaa 
-'4 



jaa 

-'5 



rpaa 
-^2 







(AA) 


(27r)'^ 


g2 (g2 _ 2q ■ pf,) 






(45) 


(27r)'^ 


g2 (^g2 _ 2q ■ p) 




r\q^) 


(46) 


{2tiY 


g2 - 2g ■ p£ 


d'^q 




(47) 




— 2g ■ p 


d'^q 




(48) 


{2nY 




d'^q 
{2nY 


r^iq^) 


(49) 


d'^q 




(50) 


{27lY 


g2 - 2g ■ p£ 


d'^q 




(51) 


{2tiY 


g2 (^g2 _ 2q . pf^ (^g2 _ 2q . p^ 


d'^q 


Jaa (^^2) 


(52) 


{2tiY 


(g2 - 2g ■ Pi) (g2 - 2q ■ p) 



The evaluation of these integrals can be done analytically and is reported in Appendix \M 



4.4 K decays 

The procedure outlined above remains true for the analysis of K iu. In the convolution 
kernel one has to simply replace p — > px {p^ = rn^ — pj^ = mj^). The form factors 
Vi and Ai remain unchanged, while in A2 one has to replace Fy^{q'^) — > Fy^{q^) = 

1 + 2HxK{(f ) + -f^7r7r(?^); whlch again amounts to the interchange ^ m^. As a 

2 4 

consequence, the full result for ^ [K Iv) can be obtained from the pion case by 
interchanging everywhere m^r with nix- 
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5 Virtual-photon corrections: results 



We collect here the results for r« 



in 



Since R, 



oc r" 



r^^P^" , we system- 



atically neglect m^-independent contributions to r| ^ that would drop in the difference. 



We also introduce the notation 

2 



me 



niK 



rriK 



(53) 



2 ? 



5.1 Leading order: 

The one loop virtual photon contributions read [TT 

'l + zi 



a 



+ 



a 



7, < , 2 
- log — + log - 

2 fi^ 1 — Zi 



log; 



1 1 



21- z, 



ze .2 



(54) 



Note that the dependence on the renormalization scale fi drops in Re/fi- Moreover, the 
dependence on the IR regulator disappears once the effect of real photon emission is 
included. 



2 4 

5.2 Next to leading order: r| ^ 

Using the notation Lg^io = (47r)^ Ll^Q^fi) and 

following expressions for the divergent and finite parts of the 0{e^p^) amplitudes 

.2 



log ^ with a = vr, K, we find the 



Vi 



(47r)4 w 9 F2 ^ 47r ' 
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(47r)4 
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F2 



47r (47rF)= 



13 - 16 

— + 8 Lg + — 

9 ^ 9 



2'^ 
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13 7 

h 

9 3 



- log(2f V^) + fi{ze) + - fi{zi) 



(55) 
(56) 



(57) 
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(47r)4 
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1 ( 5 



w V 12 



4Lc 



F2 



+ 



+ 



4^ (47rF)2 



19 ^ .r 
+ 4L9 

36 V 2 ^ 



DO 3 
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(5J 



In terms of the building block functions En{x), Rn{x), T'^'^(x), and T^^{x,y) defined in 
Appendix [XI (Eqs. 11291 and 1149111511) . the finite functions /i,2,3 read: 



/i(^) 



f2{x) 



f3{x,y) 



97 4 / / ~ 1 

— + - ( 4 iRoix) + - logx ) - Riix) - AR2{x) + R^{x] 



+ - [-9Eq{x) + 6Ei(x) + SE2{x) - GE^ix) + E4(x; 

490 + 3(147 - 327r2) 1_, , 1 T^^(x) - T^^fO) 

= ^ + -T'^^ix) + — — + R2{x) - Ro(x) 

108 3 ^ ^ 3 X iv ; ov ; 



+ 



(59) 



1 — X 



X ( R2{x) - Ro{x) ) + 2^ ( -^2(3;) - Eo{Xj 



+ X 



540 + 3(147 - 327r2 
108 



- ( 2(Eo(x) - E2ix)) + Esix) - Ei{x] 



(60) 



25 , ^r^K.. I l T^^(x,i/)-T^^(0,y) 
" 108 + 6 ^^'^^ + 6 x/y 

+ ^ (4 - y) (i?2(x) - /2o(x)) + [my) - Eo{y) 



1 1 

+ 77- 



61 — x/y 



E^{x) ~ E^{x) - E^{y) + EM 



+ x/y (4 - y) i?2(x) - Ro{x) + (x/y - 2) E2{x) - E^{x) 



(61) 



Note that the functions fi^2{x) and f^{x^y) are non-singular for x ^ (corresponding to 
rrii 0). 

In order to make the O(e^p^) amplitude UV finite we introduce in the EFT a local 
counterterm. By power counting such a term cannot distinguish K and tt decays. Its 
contribution to the amplitude is: 



2 m\ {iicf"" 
CT~^ F2 (47r)4 



d2 ^ df^+df\lJi) ■ 

— + + rcTW 



w 



w 



(62) 
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with 



d2 = (63) 
d^' = (64) 
df\li) = -hM + fhoif^). (65) 
The finite couphng rcrif^) satisfies the following renormalization group equation: 

/x|^rcT(/i) = - (4<) + 2 4^)(/i)) . (66) 

6 Matching 
6.1 Strategy 

2 4 

Within ChPT, the loop calculation of ^ produces an ultraviolet divergence proportional 
to (a/7r)m^/(47rF)^, indicating the need to introduce a local operator of O(e^p^), with 
an associated low-energy coupling. While the divergent part of the effective coupling is 
fully determined by our loop calculation, in order to estimate its finite part one needs to 
go beyond the low-energy effective theory and use information on the underlying QCD 
dynamics. 

In full generality, the 0{a) virtual-photon correction to the tt ^ £u amplitude is given 
by a sum of contributions that share the following convolution structure: 



2 



QCD 



(67) 



where is a known kernel, Hqcd stands for one of the invariant form factors appearing 
in Eq. [2H and one has to expand the r.h.s. up to O(e^p^) in the chiral power counting. 
In the framework of the low-energy effective theory, when calculating T'^ ^ we use the 

O(p^) ChPT representation for the form factors, IIqcd ^^ckpt EZl While this 

representation is valid at scales below (and generates the correct single- and double-logs 
upon integration) it leads to the incorrect UV behavior of the integrand in [671 which is 
dictated by the Operator Product Expansion (OPE) for the iyVP) and {VAP) correlators. 
As anticipated, this forces the introduction of a local operator of 0{e^p'^) whose finite 
coupling is a priori unknown, so that: 



rjie 



ChPT 



2 4 2 4 

The physical matching condition ^ \chPT = ^ \qcd in principle allows one to deter- 
mine the finite part of the counterterm. From the above discussion it is evident that the 
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counterterm arises from the UV region in the convolution of Eq. [671 so in order to estimate 
it we need a suitable representation of the correlators which is valid for momenta beyond 
the chiral regime. This poses a complex non-perturbative problem that we are not able to 
solve within full QCD. 

The problem becomes tractable if we work within the context of a truncated version 
of large-iVc- QCD, in which we replace Uqcd ^qcDo^ and Uchpt ^ ^chPT^- In this 
framework we approximate the full QCD correlators by meromorphic functions, i.e. we 
assume that the correlators are saturated by the exchange of a finite number of narrow 
resonances (at large Nc one would have an infinite number of resonances). The rele- 
vant resonance couplings are fixed by requiring that the correlators obey suitable sets of 
QCD short-distance constraints [23] (see discussion in the next section for an assessment 
of the model-dependence). Correspondingly, in the chiral effective theory the correlators 
are obtained by considering only tree-level diagrams involving Goldstone modes, with the 
couplings of higher order operators (in our case Lg and Lio) consistently determined by in- 
tegrating out the resonance fields. In this framework we are able to perform all integrations 
analytically and we determine the local coupling by the matching condition: 

Note that since we are using the large-iVc representations for the QCD and ChPT form 

4 

factors (IIqcDoc. and H^^p^^^), our matching procedure is going to miss corrections to 
c^^(/i) sub-leading in the l/A^^c expansion, which are responsible for the "double- log" 
scale dependence of the counterterm. 



6.2 Meromorphic approximation for the form factors 

In order to implement the program described above, we need a suitable representation of 
the hadronic correlator of Eq. [2ll to be used in the convolution integral of Eq. [22l In 
Refs. [211 [251 [261 [271 [28] one can find analysis of the tyVP) and AP) Green Functions, 
describing them in terms of simple meromorphic functions that respect the constraints 
imposed at low-momentum transfer by chiral symmetry and at high momentum-transfer 
by the OPE. For correlators that are order parameters of spontaneous chiral symmetry 
breaking, such as I^VP) and (^AP), this is a sensible approximation well supported by 
a number of studies [23] . 

Using the LSZ reduction formula, it is simple to extract the form factors from the 
correlators of Refs. [261 [27] . Denoting My and Ma the masses of vector and axial- vector 
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meson resonances, and using W = p — q, we find: 



^^^^'^^ - {q^ - M'y)iW^ - Ml) ^^^^ 
"^^^^'^^ - (g2-M^)(1^2-Ml) ^^^^ 

To leading order in powers of and W^^, tlie above results reproduce tlie CfiPT results 
to 0{p^), Eqs. [271 Ell Ell [301 provided one identifies the low-energy constants with their 
resonance-saturated values Lg F^/(2My) and Lio -F^/4(1/My + 1/M|), and pro- 
vided one neglects the chiral loops. Note that ^43(5^, W"^) is not relevant for our matching 
procedure, since it starts to contribute to our amplitude to O(e^p^). The dimensionless 
constants 62,3 and d2 [27] are a priori unknnown and can be fixed by imposing constraints 
on the asymptotic behavior of the Green Functions. Different results exist in the literature, 
corresponding to different choices of the resonance content of the meromorphic ansatz and 
consequently different sets of QCD short-distance constraints. These different choices will 
allow us to quantify at least in part the mo del- dependence of the final answer. Let us 
briefly discuss the two choices: 

1. The authors of Refs. [25l [26] include in their hadronic ansatz for the (VAP) cor- 
relator only the lowest lying V and A resonances and after imposing short- distance 
constraints they find 

b2 = l, h = -\, d2 = -l. (75) 

2. On the other hand, the authors of Ref. [27] include also one multiplet of pseudoscalar 
(P) resonances in the truncated spectrum. After imposing a larger set of short 
distance constraints they find 62 = 1? ^3 = 0, ^2 = 0. 

62 = 1 , &3 = , ^2 = . (76) 

For the present application it is crucial to check that the vertex functions {rvv)'^uiq^P) ~ 
Jd^x{0\T{V^ix)V!;iO)\n-{p)) and {TvA)fMp) = J d''x{0\T{V^ix) Ali0)\7r%p)) satisfy 
the correct asymptotic behavior dictated by QCD for g ^ 00. With our ansatz Tyv 
satisfies the leading and next-to-leading power behavior (0(g~^) and 0(g~^)) required by 
QCD. Concerning Tav: the ansatz of Ref. [23, [2S] reproduces the QCD behavior to 0{q~^) 
and 0(g~^). On the other hand, the ansatz of Ref. [27] has the correct 0{q^^) behavior 
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but to 0(g~^) gives a result that is twice the QCD one. Due to these considerations, in 
our analysis we will use choice 1. above for the (VAP) form factors and use difference in 
the results from choice 2. as an indicator of the model dependence. We will find that the 
spread in result is minimal, showing that the convolution integral is dominated by low and 
intermediate virtualities. This feature is quite welcome in that it makes our results more 
robust. 



6.3 Results 



The matching calculation is straightforward but tedious. It involves (i) inserting the large- 
Nc form factors of in the convolution representation of Eq. [221 (h) reducing the resulting 
integrals to scalar Passarino-Veltman functions; (iii) expanding the full result in powers 
of m^^^/My, up to order (m/My)^; (iv) finally, subtracting the ChPToo result from the 
expanded full result, thus obtaining the counterterm amplitude according to Eq. [691 The 
details of this calculation are reported in Appendix [Bl 

Using the coefficients 62,3 and ^2 as determined in Ref. [26], and defining za = Ma/ My, 
we find 



+ 



a 



184 



-1 + 4)^ 



[-37 



7 

34 



11 

y 



log 



3l4 - 17zl 



Ml 
114] 



19, 



34 (-1 + z 



2\i 



[-7 + 54 + 4-4] log^A 



(77) 



If one uses instead the values of 62,3 and ^2 from Ref. 
obtained by adding to Eq. [77] the following expression: 



the counterterm amplitude is 



2 a ni f 



log 2: 



47rM2 3(-l +4; 



(7J 



We defer a full discussion of the implications of this result to Section [HI Here we wish 
to point out that our matching procedure captures in full the " single-log" scale dependence 
of the counterterm as dictated by the renormalization group. This means that the scale 
dependence of Eq. [77] cancels the bulk of the scale dependence from chiral loops, leading 
to a very stable result. 
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7 Real-photon corrections 



7.1 Radiative decay in ChPT 

The amplitude for the radiative decay 7r~^(j>) (^~^iPe)^{Pu)l{<l) can be written as [20] : 



rprad 

L 



i2eGFF^V:,el{q)[B^ 



H"" I, 



_2p ■ q 2pi- q 

iVie^'^^^qc.Pp -Ar{q-{p- q) g'^'' - {p - qf q' 
ul{,Pu)1vv{pi>) , 



(79) 

(80) 

(81) 
(82) 



with the form factors Vi and Ai given to O(p^) in Eqs. [27] and [281 The part of the 
amphtude proportional to i?'' is referred to as "Inner Bremsstrahlung" (IB) component, 
while the part proportional to H^^ is called "Structure Dependent" (SD) component. IB 
and SD components are separately gauge invariant. The radiative decay rate has a term 
coming from the IB amplitude squared, a term from the interference of IB and SD, and 
finally a term proportional to the SD amplitude squared. To the order we work in the 
chiral expansion, only the first two terms have to be considered in principle, and lead, 
respectively, to 5^ ^ and (5| ^ in the expression for i?e/^ in Eq. [91 
Introducing the dimensionless kinematical variables 



X 



2p ■ q 



'^P-Pi 



mt 



mt 



^3) 



the differential radiative decay rate is [20 
rf2r(7r^^z/7) arW(7r^£z/; 



dx dy 
fiB{x,y) 

//ATT (^5 y) 



2tx 



(1 - z,Y 
l-y + Zi 



fiB{x, y) + m^^{ Vi fl^^i^, y) + Ai fj^'^ix, y) 



{A) 



x'^{x + y — 1 — Zi) 
x{l - y + Zi) 



+ 2(1 -x)(l - ze) 



2xze{l - zi>) 
X + y — 1 — z, 



X + y 
1 (1 



1- Zi 

y + ze) 



[2zi - x"^ + 2{1 - x){l - X - y 



X X + y — 1 - 

The total rates are obtained by integrating over the physical region 



1 - X + 



(84) 

(85) 
(86) 
(87) 



2y/z^ < X < 1 — Zg + z^ 
< y <l + Zi . 



1 — X 



When integrating the IB component over the whole physical region, an infrared divergence 
arises. It must be regulated in the same way as in the virtual corrections (in our choice 
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by giving an infinitesimal mass to the photon), and it will eventually disappear when one 
calculates the observable inclusive rate. All integrals can be done analytically, and the 
results are reported in the next section. 



7.2 Results for 6^^,61^, and ^ 

2 2n 

The IB contribution to the radiative rate reads [IHj (recall the definition of 5^ ^ in Eq. [T]): 

a I zi{10 — 7zi 



eV 



71 



4(1 - zey 



15-21ze 1 
8(1 - zA 1 



Li2{l- Ze) 



where 



+ 



2 + ilLii log zi 



1- Ze 



1 3 

log ^ - log(l - Ze) - - log ze + - 



Li2{x) = - / - log(l-t) . 
Jo ^ 



(89) 



(90) 



The above formula refers to the fully photon-inclusive radiative rate. If one considers only 
the radiation of soft photons with E^^^ < u <^ one finds, up to terms suppressed by 



a \ 1 + Ze 



^^SZe - -7z r log Ze - Li2{l - ze) 

4(1 - Ze) 1 - Ze 



+ 



1 + , 
2 + ^log; 



1 "^7 



(91) 



The interference between IB and SD amplitude (parameterized in terms of the form 
factors V\ and A\) reads: 



27T (1 - ze)^ 
7 
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log Zi - Ze log Ze 



2 I 2 



9 2^^ + ^^+ 9 



■ log Ze - z\ log Ze 



(92) 



Classifying the various terms according to their behavior with the lepton mass, one obtains: 

,2 



5f 



all. 17^^ \ 2 ^ I ^ \ , 

2? U-^'-ls''' '"- + S<^'-'^'*1^'°^^' 



a 

+ — m 
27r 



4(l-2^) + (9-5z^)log2;^ 



(93) 
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Finally, we report here the purely SD contribution to the radiative rate, which is down 
by one order in the chiral expansion but does not suffer from helicity suppression. We find: 



a 
Stt 



30 



11 
60 



20(1 - 



(12 - 3z£ 



20 z^ log zi) 



(94) 



8 Phenomenology of Re/^ 



We now put together all the results obtained so far. The starting point of our phenomeno- 
logical analysis of -Rg/^^"* is Eq. |8l which organizes the electroweak corrections to the leading 
order result of Eq. [3] according to the chiral expansion. Incorporating the effects of leading 
higher order logs [7] of the form a" log"" (m^/mg) through the correction A^i, we can write 
Eq. E] as: 

(95) 



^(0),(P) 



1 + A^fi + A\^> + A 



(P) 



g2p2 "T ^g2p4 "T ^g2p6 "I" 



1 + A 



LL 



The leading electromagnetic correction in ChPT corresponds to the point-like approx- 
imation for pion and kaon [TJ [TTl [19] : 



a(P) 



F{z) 



a 
vr 



2 ^ 



m7 



mi 



- 2- 



1 + z 



13 - 19z 



1 



8(1 -z) 
Li2{l~z) 



8-5^ 
4(1-^)2 



2;logz 



(96) 



l + z 
l-z 



\ogz log(l - z) 



(97) 



The structure dependent effects are all contained in Ae2p4 and higher order terms, which 
are the main subject of this work. Neglecting terms of order (mg/mp)^, the most general 
parameterization of the NLO ChPT contribution can be written in the form 



AT^4 



log- 



m 



m 



a m 



P ziP) 



c 



IT mt 



log- 



m 



(9J 



which highlights the dependence on lepton masses. The dimensionless constants 3 do 
not depend on the lepton mass but depend logarithmically on hadronic masses, while 



c^^\m^/mp) — > as — s> 0. (Note that our c^2l do not coincide with C2,3 of Ref. [7], 

because their C3 is not constrained to be m^-independent.) 

fp\ (p\ 

Finally, let us note that the results for 3 4 and depend on the definition of the in- 
clusive rate T{P — ^ Ivi^A). The radiative amplitude is the sum of the inner bremsstrahlung 
{Tib) component of 0{ep) and a structure dependent {Tsd) component of O(ep^) [20]. The 
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experimental definition of -Rg/^j is fully inclusive on the radiative mode, so that A^2p4 re- 
ceives a contribution from the interference of Tjb and Tsd- Moreover, in this case one also 
has to include the effect of A^SpB oc \Tsd\'^, that is formally of 0{e'^p^), but is not helicity 



suppressed and behaves as Ae2r 



On the other hand, the 



usual experimental definition of R^^l^ is not fully inclusive on the radiative mode. It corre- 
sponds to including the effect of Tib in A^^a (dominated by soft photons) and excluding 



altogether the effect oiTso'- consequently c^T^ 7^ d 



8.1 Results for i?^'!^ 

Recalling the definitions Lg = (Air)'^ Ll{fi) , ip = log(mp//i2) (^ is the chiral renormalization 
scale), 7 = y4i(0, 0)/Vi(0, 0), ze = (m^/m^)^, we find: 







(99) 



--7 + 4L9+f--2Z9 + -Ek] 4 + + -iK + -4 
24 ' ^ \36 12 12 ^ 18 8 ^ 



+ 



I - log 4 + (2 + 2 - I7) log 4 + Ki-^ (0) 



+ C3^^(/i) (100) 



cf^(m^) 



4(1 - z^) + (9 - 5z^)logz^ ) + 27 1 - 2^ + ^^og 



(47rF)2 [3(1-2^)2 L 

where k^'^^ is related to the 0{p^) pion charge radius by: 

1 „ 1 „ 1 



(101) 



=4^9 tA- 

6 3 2 



In the above equations we have used the definition 



1 ^ ^ ^ 8 m't 4 

fii^i) + 2/i(^^) + /2(^f) + fsize, z^) - - lor 



--0 ~ iog — — 

9 ° mi 9 m2 



is: J 



(102) 

. (103) 
0.025 and 



The function K^'^\mi) does not contain any large logarithms {K^'^\m^) = 
i^W(O) = 0.085) and gives a small fractional contribution to cil. 

Full numerical values of c^2 3 4 ^^2^ reported in Table [H with uncertainties due 
to matching procedure and input parameters (L9 and 7 [30]). We now discuss the results 
obtained and make contact with the previous literature. 

• We find (i^^ = in accordance to a theorem by Marciano and Sirlin [21] . This result 



arises from an exact cancellation of virtual photon contributions proportional to Vi 
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iP = n) 


{P = K) 


^2 
^3 




5.2 ±0.4i3 ±0.01^ 
-10.5 ±2.3i^ ±0.53^3 
1.69±0.07Lg 


(7.84 ±0.07^) X 10-2 
4.3 ±0.4Lg ±0.01^ 
-4.73 ± 2.3m ±0.28^3 
0.22±0.0lLg 



Table 1: Numerical values for Cn of Eq. [9H1 for P = ii, K. The uncertainties correspond 
to the input values = nip) = (6.9 ± 0.7) x 10-^ 7 = 0.465 ± 0.005 [30], and to the 
matching procedure (m), affecting only c^^\ 



and Ai and real photon contribution due to the interference of SD and IB amplitudes. 
The cancellation occurs only when the fully inclusive rate is considered. 

The coefficient c^^ is a parameter-free prediction of ChPT to this order. It involves 
the O(p^) LECs Lg and Lio, related to the pion charge radius and the ratio of axial- 
to-vector form factors 7 measurable in the radiative pion and kaon decay. 

The coefficient c^r^^ receives a predictable contribution from loops in the ChPT frame- 
work, as well as a local contribution that cannot be predicted in the purely EFT ap- 
proach, denoted by c^'^{n). Both contributions are renormalization-scale dependent, 
while the sum is not. c^'^{fi) is related to the low energy coupling Tct^Ij) introduced 
in Eq. [62] by rcrifi') = — 2(47rF)2/m^ c^^{fi). Our matching procedure gives for the 
counterterm [za = {Ma/ My) and taking My = nip): 

CTr , 19m2 / Ami 7 + ll4\ ml 

+ 



9(47rF)2 V3(4vrF)2 5^2 J ^ ^2 

6 
A 



37 -314 + 174- 11^6 



364(1-4)2 

- '^1^:1^3- log... (104) 

Numerically, using za = a/2 [32], we find c^'^{mp) = —1.61, implying that the coun- 
terterm induces a sub-leading correction to C3 (see Tabled]). The model dependence 
due to different choices of the hadronic ansatz (Ref. [26] vs Ref. [27]) is negligible, 
being Ac^'^ = 0.12. The scale dependence of c^'^{fi) partially cancels the scale 
dependence of the chiral loops (our procedure captures all the "single-log" scale de- 
pendence). Taking a very conservative attitude we assign to C3 an uncertainty equal 
to 100% of the local contribution (|Ac3| ~ 1.6) plus the effect of residual renormal- 
ization scale dependence, obtained by varying the scale /x in the range 0.5 1 GeV 
(IAC3I ~ 0.7), leading to a4^'^^ = ±2.3. 
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• Finally, the coefficient c\ can be calculated in terms of the LEG Lg and the lepton 
and meson masses and decay constants. Not surprisingly we find that this effect is 
only marginally important. 

As a check on our calculation, we have verified that if we neglect c^^ and pure two-loop 
effects, and if we use Lg = F^/(2My) (vector meson dominance), our results for 02^^34 are 
fully consistent with previous analyses of the leading structure dependent corrections based 
on current algebra [3 E]- Moreover, our numerical value of A^2p4 reported in Table [2] is 

very close to the corresponding result in Ref. [7j, namely Ag2p4 = (0.054 ±0.044) x 10~^ [7] 

versus A^2p4 = (0.053 ± 0.011) x 10^^ (this work). Therefore, as far as -Rg/^ is concerned, 
the net effect of our calculation is a reduction of the uncertainty by a factor of four. 

8.2 Results for R^^P 
In the case of K decays we find: 
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where {r'^)^^ is the O(p^) kaon charge radius and 
Moreover the function K^^\mi) is given by: 
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(107) 
(108) 

(109) 

(110) 
'm,,) 



As in the pion case, the function K^^\m£) does not contain any large logarithms (K^^^ 
0.93 and irW(O) = 1.05) and gives a small fractional contribution to cv^^. 

Note that apart from missing contributions from the SD radiation, the 4 and are 
obtained from the c^l 4 and c^"^ by interchanging m^r with rriK everywhere (the underlying 
reason is given in Sect. 14.41) . The numerical values of 4 and c^^^ are reported in Table [H 
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(P = 7r) 


(P = K) 


Ag. (%) 

Aga (%) 

(%) 


-3.929 
0.053 ±0.011 
0.073 
0.055 


-3.786 
0.135 ±0.011 

0.055 



Table 2: Numerical summary of various electroweak corrections to R^^f^^ ■ The uncertainty 
in Ae2p4 corresponds to the matching procedure. 

8.3 Resumming long distance logarithms 

At the level of uncertainty considered, one needs to include higher order long distance 
corrections [7], generalizing the leading contribution Ae2p2 ~ —3a/7r log m^/mg ~ —3.7%. 
The leading logarithms can be summed via the renormalization group and their effect 
amounts to multiplying P^^j by 1 ± A^^, with [7] 



A 



LL 
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3 TT 



log- 
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3a 
7r 



log 



me 



1.00055 . 
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8.4 Discussion 

In Table [2] we summarize the various electroweak corrections to Pg/^^"*- Applying these we 
arrive to our final results: 



P 
P 



(^) 

e/jj. 
ejii 



(1.2352 ±0.0001) X 10~^ 
(2.477 ±0.001) X 10"^ . 



;ii2) 
;ii3) 



The uncertainty we quote for Pg/^ is entirely induced by our matching procedure. However, 

in the case of -Rg/^^ we have inflated the nominal uncertainty arising from matching by a 
factor of four, to account for higher order chiral corrections, that are expected to scale as 
Ag2p4 X m\l{^TiFf. 

Our results have to be compared with the ones of Refs. [7] and [S] , which we report in 
Table O While P^^^ is in good agreement with both previous results, there is a discrepancy 

in Pg^^ that goes well outside the estimated theoretical uncertainties. We have traced back 
this difference to two problematic aspects of Ref. [S]. (i) The leading log correction A/,/, is 
included with the wrong sign: this accounts for half of the discrepancy, (ii) The remaining 
effect is due to the difference in the NLO virtual correction, for which Finkemeier finds 
) = 0.058%. We have serious doubts on the reliability of this number because the 
hadronic form factors modeled in Ref. [8] do not satisfy the correct QCD short-distance 
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10^ ■ i?^"^ 


10^ ■ ' 


This work 
Ref. [7] 
Ref. [8] 


1.2352 ±0.0001 
1.2352 ±0.0005 
1.2354 ±0.0002 


2.477 ±0.001 
2.472 ± 0.001 



Table 3: Comparison of our result with the most recent predictions of R^^f^^ ■ 

behavior. At high momentum they fall off faster than the QCD requirement, thus leading 
to a smaller value of Ag^4 compared to our work. 

9 The individual niK) tv^ modes 



The approach followed in this work is designed to obtain the ratio of -niK^ — > ev^. and 
7r(i^') — s> [iVy^ decay rates, because we have neglected all the Feynman diagrams in which 
the photon does not connect to the charged lepton. Including these diagrams in ChPT 
would generate new finite parts and UV divergences, and the corresponding local couplings 
would have to be evaluated within the l/A^^c expansion described earlier. We leave this 
task for possible future work. 

However, despite the fact that we have not performed a full O(e^p^) calculation of 
T\{K) — > ivg, our results can still be used to update the theoretical analysis of these 
individual decay modes. Here we closely follow the analysis of Ref. [7]. Including all 
known short- and long-distance electroweak corrections, and parameterizing the hadronic 
effects in terms of a few dimensionless coefficients, the inclusive P — iviiYi\ decay rate 



Fp^,, , can be written as: 



Fp,,, = F(°) X <! 1 ± ^ log ^ U ^ + - F{m\lnil) 



vr rrtp \ \ tt 









1 " 







3 VTLp 

o log 

2 VTip 



+ cr ± 4 log -2 + 4'^ + cT\m,/m,)] - 4 i'' log ^ 



;ii4) 



where F*^°^ is the rate in absence of radiative corrections (see Eq. [5]), the first bracketed 
term is the universal short distance electroweak correction, the second bracketed term is 
the universal long distance correction (point-like meson), and the third bracketed term 
parameterizes the effects of hadronic structure. The function F{z) and the constants €"^34^ 

(and £2^'') already appear in R^^^ and their expressions and numerical values have been 
reported in the previous section. The only additional ingredient needed to predict the 
individual rates Fp^^i^] is the structure-dependent coefficient c\ , which does not depend 
on the lepton mass and starts at O(e^p^) in ChPT. The explicit form (for both P = vr, K) 
is given in [11] (Eqs. 5.11 and 5.14) and it depends on a combination of EM LECs of 
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O(e^p^). These have been recently estimated in Ref. [33] in the same large- A'^c' framework 
adopted here, with the final result: 



SK) 



-2.56 ±0.5 
-1.98 ±0.5 . 



;il5) 
;il6) 



So at the moment all the structure dependent coefficients Cn'^ are known to leading order 
in their expansion within the chiral effective theory (which is O(e^p^) for c^'^ and O(e^p^) 
for the other coefficients) . The 0{e^p^) contribution to cf'^ has not yet been calculated 
(this could be done by employing the techniques presented in this paper). On the basis of 
power counting we expect c\ le^pi <^ 0.5, which is consistent with the uncertainty assigned 
to c[^^ [33]. 

Finally, we discuss here a quantity of interest in the experimental analysis of Kf.2/ K^2i 
namely the K ^ iu rate with inclusion of only soft photons {u <^ ttlx)'- 



;ii7) 



Using our results on the emission of soft photons (Eq.[9T]), it is simple to show that Tk^^^^^^ (u) 
is given by Eq. 11141 provided one replaces F{z) —>■ F^°^^{z; u), with {z = m1/m\): 



-isoft , 



3 3, 2z 
z-u) = -- + -logz-^ 



1 ± 2 , 

2 ± log z 



log 



log;z 
2uo 



l + z 



Li2{l-z) 



:ii8) 



niK 



10 Conclusions 

In conclusion, by performing the first ChPT calculation to O(e^p^) and a matching cal- 
culation of the relevant low energy coupling, we have improved the reliability of both the 
central value and the uncertainty of the ratios R^^/'^^- Our final result for -Rg/^ is consistent 

with the previous literature, while we find a discrepancy in which we have traced 

back to inconsistencies in the analysis of Ref. [8]. Our results provide a clean basis to 
detect or constrain non-standard physics in these modes by comparison with upcoming 
experimental measurements. 

As a byproduct of our main analysis, we also updated the expressions for the radia- 
tive corrections to the individual 7r{K) — > iui modes, which can be used to extract from 
experiment the combinations Vud and Fx Vus- 

Finally, it is worth mentioning that the ideas and techniques discussed in this article 
can be applied (i) to perform a full O(e^p^) analysis of the individual 7r{K) — > iui modes; 
(ii) to deal with other processes that involve one pseudo-scalar meson and a lepton pair, 
such as r — > Ki't-['j]. In this case chiral effective theory techniques are not adequate, but 
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the calculation based on the large-A^c representation for the (VAP) and (VVP) remains 
adequate. 
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A Two-loop integrals 
A.l Procedure 

In order to calculate the genuine two-loop integrals listed in Sect. 14.31 above, we use the d- 
dimensional dispersive representation of the function J°'"-{q'^) [211 [22], which is easily derived 
from Eq. [36l Re-expressing all dimensionful parameters in units of m^, one obtains: 

J-(g^) = (119) 



ds r (I) 



where q = q/ma and s are dimensionless variables. Upon inserting the representation of 
Eq. 11191 in the expression for one immediately sees that the calculation is naturally 
separated in two steps: (i) a one-loop diagram involving one propagator of mass s; (ii) 
integration of the result over the variable s, with measure given by 

In order to exemplify the procedure, we report here the calculation of j^^^'^'^ _ The other 
integrals can be worked out with similar techniques. We have found extremely useful the 
results of Ref. [22] • The case considered in that paper is slightly easier, because they only 
have one mass scale in the loops (m^), while we have two. 

Inserting the representation of Eq. 11191 in the definition of j^^^'^'^ ^ and re-expressing all 
momentum variables in units of rriT^, one arrives at (recall d = 4 + 2w): 

Here q and pe are dimensionless momentum variables (to avoid clutter we are not using 
the q,pi notation) and = {rrii/m^y. Gombining the two denominators with the usual 
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trick one gets: 



jW- ^ _^4«, I nil I dx 



[ds] /"^ , f d'^q 1 



s 



7 (2vr)'^ [(g_a;p,)2_^(3;^^)] 



2 



s 

where zix^ s) = zex"^ + s(l — x) and [22 



-im«" / / dxFMi,s)] (122) 



Explicitly one has 

F^[z] = C{w) ^-+2- ^^"^ , n > 1 (124) 

with C{w) = l/(47r)^+'". Most of the non-trivial integrals that we need to calculate have 
the structure of Eq. 11221 In order to make progress one needs to identify in Eq. 11221 the 
finite part and divergent part. This is accomplished by using a set of recursion relations 
that are the subject of next subsection. 



A. 2 Recursion relations 

In close analogy with Ref. [22] one can define (for m, n integers) : 

E{m,n]Zi)= /" M /" dx{l-xrFr\z{x,s)] (125) 

J A ^ Jo 

where z{x, s) = Ze + s(l — x). If Zi 1 then E{m, n; zi) — i?(m, n) defined in Appendix 
C of Ref. [22]. By use of integration by parts in the variable x and recalling the explicit 
form of Fn[z\ (Eq. 11241) . one can derive the following useful recursion relation: 

Y{n — 2 — w) 

{3 + w + m - n) E{m,n; Zi) = —Q{w + l-n) 

1 (n)l {—w) 

- nzi (^E{m,n + l;zi) - E{m + 2,n+l]Zi)j (126) 



with 



[ds] s' 



= C'iw)Ti-w)Ti-l-w-a)^l-^ (127) 

1 [—2a) 

Reassuringly, by setting Z£ = 1 one recovers the result of Gasser-Sainio [22] . 
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Eq. 11261 is useful because it allows one to express the integrals i?(m, n < 2; zi) in terms 
of known divergent quantities {Q{a)) and the convergent integrals E{m,3 : zi). Finally, 
let us provide an integral representation for E{m, 3; Ze) (obtained by setting d = 4): 



ZeE{'m,3;Zi) 



a{x) 



2(47r 



{1-xy 



Ia S s) ZiX"^ + s{l - x) 

(1 -x)™ r aix) , f aix) - 1 
^ ^ 1 + loe ^4 



dx 



x^ 



a{x) + 1 



^ 4(1 -x) 



ZtX^ 



1/2 



(128) 



(129) 
(130) 



We have checked that the integrals above are indeed convergent, although we could not 
find an analytic expression for z^^ 1. 



A. 3 Results 

We are now ready to present results for the integrals appearing in T, 



e p 



-'l 

T(7r)7r7r 
-'1 



^(7r)7r7r 



Tirn 
-'3 



-'4 



^5 



(47r)4 
i 

(4^ 



ml"" T{-w)T{-2w)T{l-w) 
(4^0^ (l + w)r(2-2w) 

T{-w)T{-2w)T{l-w) 



{An] 



2w 



(1 + u;)r(2 - 2w;) 



-TT 



- 7 



«m™ \C{w)Y{-w)\ 



Y[-\-w)Y[-\-2w) 



2z, 



V{-w)V{-2w) 
T{-2w)T{l-w) 



(1 + w)(2 + w) r(-'«;)r(2 - 2w) 



+ i 



2m2 



E^{z,) - ]^E^{z,) - E^ize) + ^M^i) 



(47r)4 

tml+^^ [C{w)T{-w)f 



r(-i 



w 



T{-w) 



3 17 

w + 

2 4 



59 



.2+4«, rr.,.„^T., .„m2 r(-i-.^)r(-i-2^) 



zmf^'" [C(tf;)r(-ii;)]^ 



^4+4ui 



[C{w)T{-w)X 



T{-w)T{-2w) 
V{-l-w) 



4+4io 



^(0,1;^^) 
. 4 + 2w ^ ' ' 



- E(0, 2; z,) - 2E(1, 2; z,) + E(2, 2; z,) 



(131) 
(132) 



(133) 
(134) 
(135) 
(136) 

(137) 
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rpniT 

^2 



1/ml 



(4vr) 



ds 



4\ 



1/2 



dx — log 



+ s(l 
x2 + s(l- 



x] 



X] 



(^)7r7r 



1 + 



(4vr) 



(is 



1/2 



dx[l + - — %] loj 

X X'^ 



X^ Zi + s(l — x) 
+ s(l — x) 



-frizi) 



(138) 



(4vr) 



-friz,) (139) 



Note that t5^2(s) can be expressed in terms of elementary functions and Spence functions. 
The full expressions, however, are not particularly enlightening. Since t^\{s) are not sin- 
gular for Zi — 0, numerical integration is stable and sufficient for our purposes. 
For the two-loop integrals involving J^^{q^) we find: 



{£)KK 



{tt)KK 



r(l)KK 



^{■k)KK 



jKK 
-'3 



tKK 



tKK 



KK 



Ti 



(47r)4 
i 

(4^ 



ml" V{-w)V{-2w)V{l 



w) 



(47r)2"' {l + w)V{2-2w) 



m 



4to 
K 



T{-w)T{-2w)T{l 



w] 



(47r)2"' (l + w)r(2-2w) 



Eq{z^) - E2{z^) 



tmj+^^ [C{w)r{-w)] 
2~ze 



Y[-\-w)Y{-\-2w) 

r(-w)r(-2w) 

Y{-2w)Y{\-w)' 



(1 + ry)(2 + w) r(-M;)r(2 - 2w) 



+ % 



2m\ 
(4^ 

r{e)KK\ 



Eoizi) - ^^i(5,) - E^izi) + ^Mzi) 



^2 

im 



[Ciw)Y{-w)r 



2+4«, ^^...^T.. ..m2 r(-l-i/;)r(-l-2ti;) 
is: 



r(-tf;)r(-2w;) 



zm^+^'" [C(tf;)r(-u;)]^ 



/f^+24m|+^"' 



2£ 



4 + 2w 



Y{-l-w) 

r(-«;) 

^(0, 1; 5,) 



5,' E(0, 2; 5,) - 2E{1, 2; 5,) + E(2, 2; 5,) 



1/m 



(47r)4 - 2^ 

1 1 
(ix — log 

X 



s 

„2~ 



4\ 



1/2 



(140) 
(141) 



1-- tris). 



(47r)4m2 ^1 



(142) 
(143) 
(144) 

(145) 
(146) 

T.^^(5,,5.) (147) 



X Zi + s{l — x) 
x'^Zj, + s(l — x) 
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Tr 



KK 



(l)KK 



(47r)4 



4\ V2 



t2"[S] 



{e)KK _^ ^ jKK I z z 



1 + - 



— ^vr Jo 



, /_ s s \ , / x^ze + s(l — x) 
dx[z^^ -] log 



(148) 



X x^ 



X^Zt, + S(l — X) 



The finite loop contributions can all be expressed in terms logarithms and combinations 
of En{x) and the following functions: 



RJx) 



EJx) 



X 



(149) 

(150) 
(151) 



A. 4 Standard form of two-loop integrals 

A generic two-loop contribution can be cast in the following standard form 
l/(47r)2+"'): 

j2-ioops ^ [C{w)T{-w)f m^"" x{d) 
x{d) = Xq + XiW + X2w'^ + 0{w^) , 



with m = rrij^ or m = rriK- Multiplying and dividing each contribution by (/ic) 

logc = (log47r - 7s + 1) 

and performing the expansion around d = 4, one finds: 

(/ic)^"' 



4ui 



j2— loops 



(47r)4 



+ + F + 0(w) 



w 



xi + 2 xo ( lo. 



m 



m 



X2 + 2 Xi ( log — + 1 ) + Xo 



— + 2 log— + 1 



6 



(C(^) = 

(152) 
(153) 

H], with 
(154) 

(155) 
(156) 
(157) 

(158) 



A. 5 Standard form of one-loop integrals 

The one loop diagrams with one insertion from the O(p^) effective lagrangian can be cast 
in a useful standard form as well. Denoting by L{d) the generic d-dimensional LEG, 
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one has: 



y{d) = yo + yiW + y2W^ + 0{w^) 

{ncf^ ( r 



(4vr) 



(159) 
(160) 

(161) 



where m = or m = and the constant F determines the RG running of the renor- 



mahzed couphng Ui^^i). Multiplying and dividing each contribution by (/xc) 
performing the expansion around d = A, one finds: 

^ + + F + Oiw) 



2w 



—loop 


(/ic)^"' 




(47r)4 






2 


^(1) = 


-(47r)2 


F = 


-(47r)^ 



w 



2 rr/ 



2 2 
-2/0 -^((6 + -^ 



log^ 



, and 

(162) 

(163) 
(164) 



{2 W r {fi)yo +{yo + yi)r) log --yoT flog-)' (165) 



The couplings of interest to us are Lg and Lio, whose divergent parts are determined by 

1 _ 1 



(166) 



B Matching calculation 

In this Appendix we report the details of our matching calculation. The intermediate steps 
of the calculation are: 

1. Insert the large- A^c form factors of in the convolution representation of Eq. [22j 

2. Reduce the resulting integrals to scalar Passarino-Veltman functions [31]. For these 
we follow the convention of Kniehl [35] . 

3. Expand the full result in powers of nif^T^/My, up to order {m/MyY. This involves 
expanding the scalar integrals Bq{p'^ ,m\,rn^ and Co(...) in powers of ratios of the 
internal masses. This is trivial for Bq, somewhat less trivial for Cq. We derived a 
representation of Co as a two dimensional integral (see Ref. [36]) and used that as a 
starting point for the heavy mass expansion. 

4. Subtract the ChPToo result from the expanded full result, thus obtaining the coun- 
terterm amplitude according to Eq. [691 
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B.l Reduction to Passarino Veltman functions 

We use the conventions of Ref. [35] for the Passarino- Veltman functions, namely: 



and 



with 



{Bo, B^, B^^} ml, ml) = j 



{l,q^,q^,qu} 



27]-2 ^q2 _ ^2 _|_ _|_ py _ ^2 _|_ 



{Co, C^, Cf,u} (p^ {p + mf, ml, m|) = 

d'^q {^.q^l,q^lqu] 

1^2 [■g2 _ ^2 _|_ _|_ pY _ ^2 _|_ _|_ p _|_ /.p _ ^2 _|_ ^gj 



5, 



Pa. Bi 

PliPu B21 - Qf^y B22 



Cmi/ = P^iPu C21 + k^ku C22 + {p^iku + kf,p^,) C23 - g^wC24 



(167) 



(168) 



(169) 
(170) 
(171) 
(172) 



For the reduction of vector and tensor integrals to scalar Passarino- Veltman functions we 
have used the relations {A.7), {A.8) and {A.9) of Ref. 



B.1.1 Tv, 

In the reduction of Ty^ we need the following tensor and vector integrals: 

/^d Y f 2 LT/2\ 



(175) 
(176) 
(177) 
(1,78) 



Using the above definitions the amplitude reads: 



Vn 



6(47r)2 

i 

6(47r)2 
z 

6(47r)2 



6K, + (< - m^) ( \/« - 
1 



M2 



{B2i{m], Ml, ml) - B2i{ml 0, m^)) - (1 - ft:)C'2i - KC21 



[1 — k)C23 — I^C23 
1 



M2 



{B22imj, M^, mj) - B22{m], 0, m])) + (1 - k)C2a + KC24 
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with 



2M^ 



ml - cv 



Cij 

Cii 



Mt = -6 Mt Vi 



V 4^2^2 V 

Cij{mlO,ml,0,mlM^) 
Ci,{mlO,ml,M^,mlM^] 



(179) 

(180) 

(181) 
(182) 



B.1.2 Ta, 

In the reduction of Ta^ we need the following tensor, vector, and scalar integrals: 



{2TxYq^ (g2 -2q-i 



q°'q 



d'^q Ai{q\W^ 



d'^q A^{q\W^) „ 



{2tiY q^ -2q- 
d'q A,{q\W^) 

- q = -&7r P 



(2vr) 



q^ 



=A,pe + A,p^ (184) 

(185) 

=E,p'^ + Eep1 (186) 

(187) 



Using the above definitions the amplitude reads: 



Ta^ = -ie^Tf I S'ai - K + ^,.("^^ - 



+ {d-2) 



Sa^ 
A, 
Alt 



A 



— H H 



(47r)2 
i 

(4^ 

— z 
(4^ 

—2 



2 " 2 

feiCo + 62 5o(0, m^, Ml) + 63 5o(m2, , m|) 

&2I ^2 



^1 ^ ^ ^1 

<-^12 ~ 



Ml 



Ml 



^ {B2,{mlMl,m^,) - B2i{ml0,mj)) 



M2 



&2 C21 



(4vr) 



<-^23 ~ 



Ml 



Ml 



b2 C23 



(188) 
(189) 
(190) 



(191) 
(192) 
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A. 



—I 



(4vr) 



<-^24 + 



Ml 



(47r)2 

i 

(4^ 



ho] a 



24 



feiCn - 62^0(0, ml Ml) + 63^1 ("^', Ml, m] 

hi 



Ml 



Ml 



62 i?iK,0,Mi) 



with 



61 = Ml{l-h2)-Ml{l + h; 



Cij 



Qj{mi,0,mi,Ml,mi,Mi) 
a,imlO,ml,0,mlMl) . 



B.1.3 Ta, 

In the reduction of we need the following vector and scalar integrals: 

A A2iq',W^) 

= ^A2 



(2^ ^2 



Using the above definitions the amplitude reads: 

ie 



Ta2 



(47r)2 

i 

~(4^ 



2Sa2 + (2 - d)Fe - 
(2 + d2)MlCo - ^2 ^oK', Ml, m\ 

(2 + ^2)^2(711 - d2Bi{ml Ml, mj 



j /Vf2 

A(2 + Ml, Ml) - B,{ml, 0, Ml)] 



(47r)2 M2 



B.1.4 Ti.^ 

The Fy-induced amplitude reads: 

^2 



+ 



(47r)2 
1 



m2 — rni 



Tf I (m^ + mj) Coiml 0, , Ml, ml ml) 
{m\ B^{ml, Ml, ml) - ml Bo{ml Ml, m^] 
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B.2 Expansion of the relevant three-point scalar functions 

We use the following representation for the Co function as a basis for the large mass 
expansion: 

Co{p'',k\{p + kY,mlml,ml) = / dx / dy , . „ —, (207) 



^0 



ax"^ + by'^ + cxy + dx + ey + f 



with 



a = [p + kf 

b = 

c = {p + kf+p^-k"^ 

d = ml — m\ — {p + kY 

2 2 2 

e = 7712 — — p 
f = mj . 

We then find (we give results up to the needed order): 

1 1 / \ 
Co{mlO,ml,M^,mlml) = — ^ log ^ - m.Mog — ^ + . . . (20J 

. 9 2 9 2 s ^ f , Ml\ ml + mi f .Ml 

C„(mlO,ml,0,,nlM-) = ^ ^1 + ,og + (l + 2l„g -| 

ml + mjml + mj f M? 



9M6 



1 + 3 log ^ 



V \ 



ml + mjmt + mfml + ml f ^ ^ ^ M^^ . ^ ^^gg) 



16 



2 

mf 



Co(m2,0,m^,M2,m2,M2) = ^ + ^ (^^m] + ^ml - m]\og^ 

^ 1/28 4 5 2 2 ^1 4 o 41 , ^ 

+ — s- — mm. H — 3m« log — ^ + 

V 9 ^ 36 " ^ 90 " ^ ^ ' 



Co(m,^0,m^,M2,m,^Ml) = log^ + (4, m^, m^) 



M6 



(210) 



+ - , (211) 
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with ZA = Ma/Mv. The functions f^^'^\z\,m},ml) have a simple but lengthy expression: 



J{4)(.2 ^2 ^2^ 



f'\z\M.ml) 



+ 



+ 



+ 



-^log^ + 



2; 



2(-i + 4)^ 



K4(i + 4) + ^K2-34 + 4)) 



log 



m^(l + 24) 



log 



+ 



m^(3-34 + 4) 



34( 



log 2^ 



mjml{A - 54 + 4) + "^^(1 + 44 + 4) 



m^(3 



3( 

154 + 104 



log 4 



64( 



■1 + 4P 



24 



74 



^t(i + 4) 

(-1 + zir 

24) 



64(-i + 



(212) 



(213) 



B.3 Results 

Recalling the definition za = Ma/Mv and neglecting as usual the m^-independent terms 



that drop in Re/fj,, we find 



rpCT 



CT 



rpCT 



,2 a 
,2 a 



2 log 



4 19 
9 9 ^ 



2 ' -1^1 log 



(214) 
(215) 



log 



Ml 



37 - 634 



214 + 5z\ + 12(7 - 10z\ + 44) logZA 



184( 



-1 + 4)^ 



(216) 



rpCT 
-'-A2 



^p2 a ml 



471 



^, M2 12-164 + 44 + 4(12-15^2 
- » log ■ 



54) log Za 



+ d' 



r 3(-i + 4)3 

8 - Uzl + 64 + 2(12 - 15z^ + 54) log^A 
3(-l+4)3 



Using the input from Ref. [27] the second line of T^'^ should be replaced by: 
37 - 1484 + 1684 - 524 - 54 + 12(7 - 284 + 274 - §4) log^A 



184( 



-1 + 4)3 



(217) 



(218) 



By comparing these expressions with the ChPT ones, one can easily verify that the match- 
ing procedure captures in full the single-log renormalization scale dependence, as expected 
from the I/Nq expansion. 
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